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TIME  CORRECTION  IN  PASSIVE  RANGING: 
BREAKTHROUGH  OR  BOOTSTRAP?* 


Introduction  and  Summary 

The  concept  of  time  correction  in  passive  ranging  is  described 
in  reference  (a).  It  was  developed  by  Daniel  H.  Wagner, 
Associates,  who  believe  it  to  be  a  breakthrough  in  the  subject 
of  passive  ranging.  The  concept  has  been  applied  to  the  four- 
bearings  problem  for  which  it  was  conceived  and  to  the  Ekelund 
range  estimate  involving  two  bearing  rates.  A  time-corrected 
Ekelund  range  estimator  is  being  considered  for  part  of  the 
target  motion  analysis  in  a  new  fire-control  system. 

Does  the  concept  of  time  correction  represent  a  significant 
breakthrough,  as  claimed,  or  is  the  apparent  gain  an  illusion? 

We  have  examined  this  question  carefully  and  have  found  no 
support  for  the  claims  that  have  been  made  for  the  method.  We 
have  concluded  that  the  method  has  no  special  value  in  passive 
ranging  and  should  be  discarded.  Our  arguments  are  presented 
below. 

The  Time-Correction  Method 

Time  correction,  as  developed  in  reference  (a),  is  a  method  for 
the  computation  of  a  range  R*  at  a  time  t*  from  the  four  bearings 
Bl’  B2  B3*  B4  °bserved  at  the  corresponding  times  t^,  t^,  t^,  t^ 
which  are  assumed  to  be  in  the  order  t^  <  t£  <  t^  <  t^.  We 
outline  the  derivation  of  the  equations  below. 


This  research  was  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-70-C-0322. 
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We  use  the  standard  naval  coordinate  system.  For  i  ■  1,  2,  3,  A, 
let 


t^  ■  time  at  the  l*h  observation 

-  bearing  of  the  target  from  the  tracking  submarine 

-  range  of  the  target  from  the  tracking  submarine 

uf^  ■  range  component  of  target  velocity  measured 
outward  along 

■  normal  component  of  target  velocity  measured 
normal  to  B^  in  the  direction  of  positive 
(clockwise)  angular  rotation 


Also,  let 


y  -  distance  (directed)  tracking  submarine  moves  in 

^  the  direction  B,  from  time  t.  to  time  t. 
k  k  i  j 

v*  -  distance  (directed)  tracking  submarine  moves 
^  normal  to  B,  from  time  t.  to  time  t. 

k  1  J 


It  will  be  convenient  to  use 


■  cj  -  V  Bu  ■  Bj  -  V  sij  • 8l"  V  cu  ■ 009  Bij 


Assuming  linear  target  motion  it  is  easy  to  derive  the  equations 
S, 


r2  c 


12  "  ^1  fc12  "  12 

(1) 

3 

34  "  %3  C34  "  V34 

(2) 

'24  "  R4  “  Ur4  C24  +  y24 

(3) 

C13  "  “bS  +  Url  S13 

(4) 

These  equations  are  equivalent  to  the  equations  (2-5)  (2-6)  (2-7) 
and  an  unnumbered  equation  on  page  2-18  of  reference  (a) . 
Eliminating  R2>  u^,  and  u^  from  the  above  equations  and  solving 


for  R,  we  obtain 
4 


R4  "  R4  +  al  url  +  a2  ur4 


(5) 
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where 


*4  "  (V34  1  *34  -  C13  V12  1  *12  “  y24  g12)  /  (812  "  e345  (6) 


al  "  S13  1  (012  _  034J»  S2  "  *24  612  /  (B12  “  034J  (7) 


612  "  *S12^*12*  (C13/C24)’  e34  "  S34  /  *34  (8) 


Equation  (5)  expresses  the  range  R^  in  terms  of  an  estimate 

A 

that  can  be  computed  from  observed  quantities,  and  two  range 
components  of  target  velocity. 

In  the  time-correction  method  the  two  range  components,  url  and 
u  of  target  velocity  are  eliminated  by  introducing  two  bearings, 
and  Bg,  at  times  t,.  and  tg  by  using  two  equations  of  type  (3) 
above.  They  are 


R5  C45  "  R4  +  Ur4  *45  “  y45 


R,  Ct ,  -  Rc  C1c  +  u  ,  t.c  -  \i.c 
6  16  5  15  rl  56  56 


The  general  equation  of  this  type  is 


Rj  Cjk  “  Ri  Cik  +  urk  *ij  '  pij  (11) 

Equations  (3)  (9)  (10)  are  special  cases  of  equation  (11). 

Eliminating  R^  and  R,.  from  equations  (5)  (9)  (10)  we  obtain  the 
equation 

R6  C16  "  (C15/C45)  (^4  "  ^45^  "  y56  +  a3  Url  +  a4  Ur4  (12) 


where  a^  and  a^  depend  on  t^  and  tg.  Let  t  and  t  be  the  values 

of  t^  and  tg  respectively  that  make  ■  0.  Also,  let  B 

and  B  be  the  corresponding  bearings  of  the  target  at  these 
* 

times,  and  let  R  be  the  corresponding  value  of  Rg.  Then 
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t  ■  (^2^2  ”  ^34C4^  ^  ^12  ”  ^34^  (13) 

t*  -  t  -  S13  cos  (B  -  Bj)  /  (012  -  e34)  cos  (B  -  B^)  (14) 

R  cos  (B*  -  B^  -  (R^  -  y45)  cos  (B  -  B^)  /  cos  (B  -  B^)  -  y^ 

(15) 

Equations  (13)  (14)  (15)  are  the  same  as  equations  (2-14)  (2-18) 

(2-20)  respectively  of  reference  (a)  when  our  subscript  sequence 

(1,  2,  3,  4)  Is  replaced  by  the  sequence  (1’,  1,  2',  2)  of 

reference  (a).  (It  Is  convenient  to  retain  the  subscripts  5  and 
4  1 

6  in  and  y ^  to  avoid  ambiguity.) 

'ft  ft 

The  range  R  is  called  "the  time-corrected  range"  and  t  is 
called  "the  best  range  time."  The  time-correction  method  consists 

it 

of  using  equation  (15)  to  compute  R  and  equation  (14)  to  compute 

•k 

the  corresponding  time  t  at  which  it  applies. 

Discussion 

Does  the  above  method  offer  a  significant  breakthrough  in  passive 

ranging?  Various  claims  are  made  for  the  method  in  reference  (a) . 

It  is  asserted  that  the  basis  for  the  method  is  that  "every 

ranging  maneuver  has  associated  with  it  a  best  range  time",  at 

which  time  the  error  in  the  range  component  of  velocity  "will 

have  minimum  effect  on  the  accuracy  of  the  range  solution."  And 
* 

when  t  is  in  the  future,  it  is  claimed  that  the  four  bearing 

4  4 

observations  and  the  three  distance  measurements  (y34»  ^45* 
and  y^g)  "provide  an  accurate  range  at  a  future  time  with  no 
additional  information  needed."  Are  these  claims  justified? 

The  derivation  is  deceptively  simple.  No  optimization  is  needed. 
Only  simple  algebraic  manipulations  of  well-known  equations  are 
required  to  eliminate  the  two  range  components,  u^  and  ur4,  of 
target  velocity  from  the  range  equation  (5).  It  is  surprising 
that  such  a  simple  derivation  yields  a  range  estimate  that  is 
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superior  to  other  range  estimates,  as  claimed.  And  it  is  even 
more  surprising  that  the  result  remained  undiscovered  so  long, 
since  the  derivation  involves  only  simple  mathematics. 

A 

Since  t  is  not  obtained  by  optimizing  a  payoff  function,  why 
is  it  called  "best"?  Apparently,  t  is  called  the  best  range 

A 

time  because  the  corresponding  range  R  has  been  written  in 
equation  (15)  in  a  form  that  does  not  involve  a  range  component 
u^  of  target  velocity.  Immediately  following  the  equation,  which 
is  numbered  (2-20)  in  reference  (a),  it  is  stated  that,  "The 

A 

range  R  is  the  time-corrected  range.  Note  that  all  of  the 
terms  in  Equation  (2-20)  can  be  measured  by  the  SSK.  The  time 

A 

t  at  which  this  applies  is  called  the  best  range  time... 

Can  all  of  the  terms  in  equation  (15)  above  be  "measured"? 

—  A 

Certainly,  B  and  B  can  t  be  measured  in  the  same  way  that  B^, 

B2,  B^,  B^  are  measured  by  recording  the  observed  value  at  the 
time  of  occurrence.  The  bearings  B  and  B  are  the  bearings  at 
the  times  t  and  t  .  Obviously,  if  t  and  t  fall  outside  the 

-  ft 

tracking  interval  (t^,  t^)  ,  the  corresponding  bearings  B  and  B 
can't  be  observed.  Nor  can  they  be  observed,  if  t  and  t  fall 

—  A 

inside  the  tracking  interval.  The  values  of  t  and  t  in 
equations  (13)  and  (14)  require  the  values  of  all  four  bearings. 

—  A 

Hence,  t,  is  the  earliest  time  at  which  t  and  t  can  be  computed, 

,  ^  _a 

and  it  then  is  too  late  to  observe  B  and  B  . 


Hence,  B  and  B  must  be  "estimated"  from  the  bearing-time  plot, 

—  A 

by  interpolation  when  t  and  t  fall  inside  the  tracking  interval 
and  by  extrapolation  when  they  fall  outside  the  tracking  interval. 

— !  •  A 

In  either  case  the  errors  in  estimating  B  and  B  may  be  large. 
Here,  interpolation  may  be  as  inaccurate  as  extrapolation,  since 
the  SSK  is  maneuvering  during  the  tracking  interval  and  bearing 
observations  usually  are  not  made  during  the  turns.  However, 

A 

the  resulting  errors  in  the  cosine  factors,  and  in  R  ,  in 
equation  (15)  would  be  small.  Hence,  the  fact  that  B  and  B 
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must  be  estimated  by  interpolation  or  extrapolation  is  not  a 
major  disadvantage.  The  error  is  comparable  to  that  made  in 
using  the  small-angle  approximations. 

4 

Errors  also  will  occur  in  estimating  the  two  distances  p^,.  and 

,  since  they  also  must  be  obtained  by  interpolation  or 
_  * 

extrapolation  after  t  <-t_)  and  t  ("O  have  been  computed. 

3  0  * 

These  errors  may  produce  larger  errors  in  R  than  those  produced 

—  A 

by  the  errors  in  estimating  B  and  B  . 

More  serious  questions  concerning  the  method  can  be  raised  as 
follows: 

a.  What  is  the  motivation  behind  the  procedure?  What 
objective  are  we  trying  to  achieve?  Is  it  merely  to 

A  A 

find  an  equation  for  the  range  R  at  some  time  t  in  a 
form  that  Involves  only  quantities  that  the  SSK  can 
measure  or  estimate? 

b.  Does  the  procedure  lead  to  a  unique  solution  that  has 
some  desirable  property?  If  so,  what? 

c.  After  we  find  R  at  time  t  ,  how  will  we  obtain  an 
estimate  of  the  range  R  at  time  t  when  the  time  t  at 
which  we  need  a  range  estimate  does  not  coincide  with 

A 

t  ?  Will  this  require  that  we  use  an  equation  of  type 
(11)  above  that  requires  a  range  component  u^  of  target 
velocity?  If  so,  what  have  we  gained  by  first  computing 

A 

R  ? 

If  the  objective  of  the  method  is  to  obtain  an  equation  for  range 
in  a  form  that  involves  only  quantities  that  the  SSK  can  measure 
or  estimate,  why  stop  at  one?  There  are  many  such  equations. 

What  special  properties  does  the  "solution"  in  equation  (15) 
have  that  it  should  be  selected? 

Other  "solutions"  of  the  type  displayed  above  can  be  obtained 
by  interchanging  two  subscripts.  This  possibility  seems  to  have 
been  recognized  in  some  parts  of  reference  (a) ,  and  used  in  a 
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special  case  in  chapter  4.  But  the  implications  concerning  the 
validity  of  the  general  method  and  the  uniqueness  of  the  solution 
were  not  explored. 

Some  of  the  "solutions"  obtained  in  this  way  reduce  to  the  same 
A  A 

forms  for  t  and  R  in  the  small-angle  approximations,  and  others 
do  not.  For  example,  if  we  replace  equations  (9)  and  (10)  by 
the  equations 


'5 

6  C46 


R.  C  .  +  u  t. 

4  14  rl  45 


R5  C45  +  Ur4  fc56 


and  repeat  the  steps  of  the  solution,  the  equations  for  t  and 
a 

R  are  not  identical  to  equations  (14)  and  (15)  but  reduce  to 
the  same  small-angle  forms  that  equations  (14)  and  (15)  become. 
The  equation  for  t  does  not  do  so.  On  the  other  hand,  if  we 
interchange  the  subscripts  2  and  3  in  the  original  solution, 
we  get  a  solution  for  which  the  small-angle  approximations  for 

£  ff  _ 

t  and  R  ,  as  well  as  that  for  t,  do  not  reduce  to  the  same 
small-angle  forms  obtained  in  the  original  solution.  This 
solution  is  written  in  full  belcw: 


1  "  (613  t3  "  024  t4)  /  (B13  “  624) 


(13a) 


*  _ 


t  -  t  -  S12  cos  (B  -  B1)  /  (B, ,  -  89/.)  cos  (B  -  BA)  (14a) 


13  24 


ft  ft  A  A  ^  4 

R  cos  (B  -  B^)  -  (R^  -  y^)  cos  (B  -  B^)  /  cos  (B  -  B^)  -  y,.^. 


(15a) 


where  now 


™  ^v24  ^  fc24  ”  ^12  vi3  ^  ti3  ”  ^34  ^13^  ^  (6a) 


13  K24' 


613  "  (S13  /  t13)  (C12  /  C34)’  824  =  S24  /  fc24 


(8a) 
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fg 

It  is  evident  that  we  can't  call  the  time  t  in  equation  (14) 

the  best  range  time,  since  it  is  neither  unique  nor  optimal  in 

any  demonatrated  way.  It  may  be  a  good  range  time,  but  the 

property,  if  any,  that  makes  it  better  than  an  arbitrary  time 

has  not  been  displayed  in  reference  (a)  or  elsewhere,  to  our 

knowledge.  The  property  that  led  to  its  introduction  -  that  it 

* 

permits  us  to  write  an  equation  for  the  corresponding  range  R 
in  a  form  that  involves  only  quantities  that  the  SSK  can  measure 
or  estimate  -  is  not  a  sufficient  reason  for  its  use.  Any 
arbitrary  time  has  this  property. 

Under  the  assumption  of  linear  target  motion,  on  which  the 

derivations  above  and  in  reference  (a)  are  based,  four  bearing 

observations  (made  on  a  ranging  maneuver  that  doesn't  consist 

of  a  single  linear  path)  are  sufficient  for  a  complete  TMA. 

Hence,  the  range  R  at  any  arbitrary  time  t  can  be  written  as  a 

function  of  t,  the  four  observation  pairs  (t^,  B^) ,  i  ■  1,  2,  3,  4, 

and  components  of  own-ship  motion.  It  is  not  necessary  to 
—  * 

compute  two  times  t  and  t  and  then  estimate  the  corresponding 
—  * 

bearings  B  and  B  by  interpolation  or  extrapolation  from  the 
bearing- time  plot.  The  only  bearings  involved  are  B^,  B2,  B^,  B^. 

The  TMA  requires  only  the  solution  of  four  simultaneous  linear 
equations  in  a  set  of  target  parameters,  such  as  x^,  y^,  u^,  u^, 
where  ux  and  u^  are  the  components  of  target  velocity.  The 
solution  has  been  obtained  many  times.  In  fact,  Spiess'  form 
of  the  solution  is  discussed  in  chapter  5  of  reference  (a) .  By 
means  of  the  TMA  we  can  write  an  equation  for  the  range  R  at 
an  arbitrary  time  t,  in  a  form  that  involves  only  quantities 
that  the  SSK  can  measure  or  estimate. 

The  solution  for  R  is  simple  when  t  coincides  with  one  of  the 
observation  times.  For  example,  when  t  ■  t^,  the  solution 
becomes 
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(16) 


Rz 


V14  *24  *34  S23  ~  V24  *14  *34  S13  *  V34  *14  *24  S12 
=_ ^13^24 


where ,  as  before. 


distance  SSK  moves  normal  to  B,  in  the  interval 

<V  V 

sin  (Bj  -  Bi) 


We  do  not  recommend  the  use  of  the  estimator  (16)  for  R^,.  However, 
equation  (16)  is  a  simple  equation  for  the  range  in  terms  of 
the  four  (t^,  B^)  pairs  and  three  components  of  own-ship  motion. 
Also,  it  is  exact,  that  is,  it  will  yield  the  exact  value  for 
the  range  at  time  t^  in  the  absence  of  errors  in  the  measured 
quantities.  And,  certainly,  the  SSK  can  measure  all  the 
quantities  in  equation  (16) .  Do  these  facts  make  t ,  a  "good" 

ft  ^ 

range  time?  If  not,  what  property  does  t  possess  that  makes 

ft  * 

it  better?  Will  any  gain  that  is  obtained  from  using  t  and  R 

ft 

be  lost  in  computing  R^  from  R  ,  if  t^  is  the  time  at  which  a 
range  estimate  is  needed. 

It  is  surprising  that  nothing  is  said  in  reference  (a)  about  the 

ft  ft 

problem  of  converting  from  the  estimate  R  at  time  t  to  an 
estimate  of  the  range  at  the  time  at  which  it  is  needed.  How  is 
this  conversion  to  be  made?  By  solving  for  the  TMA?  If  so, 
won't  we  finish  with  the  same  deterministic  relationships,  but 
in  a  more  circuitous  form,  with  many  more  chances  for  errors  in 
the  measurements?  Will  the  estimate  of  by  such  a  method  be 
more  accurate  than  that  obtained  from  equation  (16)? 

A  possible  way  of  avoiding  these  questions  is  to  state  that  the 

ft 

submarine  commander  can  control  the  value  of  t  by  his  choice  of 

* 

ranging  maneuvers  and  hence  can  make  t  equal  to,  or  close  to, 
the  desired  time,  say  t^.  In  fact,  this  point  is  offered  in 
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reference  (a)  as  a  major  advantage  of  the  time-correction  method. 

* 

We  contend  that  control  of  the  value  of  t  is  not  a  practical 
solution  to  this  problem.  The  commander  of  own  submarine  has 

4r 

only  partial  control  over  t  and  that  control  is  rather  tenuous. 
The  commander  of  the  target  submarine  must  cooperate  by 
maintaining  constant  course  and  speed.  Even  when  he  does  so,  the 
problem  of  the  selection  of  maneuvers  and  times  t^,  t£,  t^»  t^ 
that  will  make  t  fall  at  t^  is  not  an  easy  one  to  solve,  since 
the  value  of  t  can’t  be  computed  until  the  last  bearing 
has  been  recorded.  It  is  the  type  of  problem  that  can  be  solved 
after  the  fact  but  is  very  difficult  to  solve  in  real  time  and 
sequence.  Even  if  it  can  be  done,  there  remains  the  problem  of 
how  to  assure  linear  target  motion  in  the  relevant  time  interval. 

Criteria  For  Selection  of  a  Range  Estimator 

The  equations  derived  above  have  been  based  on  two  important 
assumptions,  as  follows: 

(1)  The  target  submarine  maintains  constant  course  and 
speed  throughout  the  relevant  time  interval,  which  is 
the  maximum  interval  spanned  by  the  times  involved  in 
the  computation. 

(2)  There  are  no  errors  in  the  measurements  of  bearings, 
including  those  obtained  by  interpolation  or  extra¬ 
polation,  or  in  the  measurement  of  own-submarine 
motion. 

Both  assumptions  are  questionable,  at  best,  and  we  should  consider 
the  effects  of  departures  from  these  assumptions. 

The  argument  so  far  has  been  based  solely  on  formula  accuracy 

* 

under  these  assumptions.  The  main  argument  for  the  use  of  R 
in  equation  (15)  seems  to  be  that  equation  (15)  is  exact  and 
contains  only  quantities  that  the  SSK  can  measure  or  estimate; 

10 
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whereas  some  other  equations  for  a  range  estimator  are  inexact, 
such  as  equation  (6)  as  an  estimator  for  R^;  or  they  contain 
quantities  that  the  SSK  can't  measure  or  estimate,  such  as 
equation  (4).  More  precisely,  equation  (4)  contains  two  factors, 
u^  and  u^,  that  the  analyst  hasn't  taken  the  trouble  to  express 
in  terms  of  quantities  that  the  SSK  can  measure  or  estimate. 

If  this  had  been  done,  equation  (5)  would  have  reduced  to 
equation  (16)  when  simplified. 

Formula  accuracy  is  only  one  characteristic  of  a  range  estimator, 
and  not  the  dominant  one  in  most  applications.  The  main 
characteristics  of  a  range  estimator  are: 

(a)  Formula  accuracy 

(b)  Ruggedness  to  target  maneuvers 

(c)  Sensitivity  to  error  measurements 

It  is  not  likely  that  an  estimator  will  dominate  all  other 
contenders  in  all  respects.  For  example,  the  estimator  for  R^ 
in  equation  (16)  is  exact,  and  hence  superior  in  formula 

A 

accuracy  to  the  estimator  R^  for  R^  in  equation  (6) ,  and  to  the 
estimator  Rg  obtained  by  dropping  the  term  in  equation  (6) , 
and  to  the  Gkelund  estimator  R^  (a  limiting  form  of  R^) ,  since 

R4  in  (16) 

probably  is  not  as  rugged  to  target  maneuvers  as  some  of  the 
other  estimators,  and  undoubtedly  is  more  sensitive  to  random 

A 

errors  than  R 

ruggedness  and  sensitivity,  or  even  the  formula  accuracy,  of  an 

* 

estimator  for  R^  that  starts  with  R  in  equation  (15) ,  since  it 

* 

depends  on  the  procedure  used  to  calculate  R^  from  R  . 

If  there  is  no  dominant  estimator,  how  do  we  select  one?  We 
believe  that  this  decision  requires  a  careful  examination  of  the 
use  that  will  be  made  of  the  range  estimate  and  the  selection  of 
a  suitable  performance  measure  or  payoff  function.  We  then 


Rg,  or  Rj,^.  It  is  impossible  to  estimate  the 


R^,  Rg,  and  R^  are  biased  estimators.  However, 
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compute  the  payoff  function  for  possible  target  maneuvers  and 
error  distributions,  select  the  dominant  estimator,  if  it  exists, 
or  use  a  weighted  mean  otherwise. 

In  many  applications  formula  accuracy  is  relatively  unimportant, 
particularly  when  likely  target  maneuvers  or  measurement  errors 
will  produce  range  errors  that  are  much  larger  than  the  inherent 
bias  in  the  formula.  Under  such  conditions  it  is  absurd  to  select 
an  estimator  on  the  criterion  of  formula  accuracy,  alone  or 
predominantly. 

Comparisons  of  range  estimators,  or  complete  TMA  estimators, 
often  (usually?)  have  been  made  on  a  partial  basis.  We  believe 
that  all  characteristics  of  the  estimator  should  be  considered. 

It  should  be  possible  to  analyze  and  compare  leading  contenders 
in  respect  to  the  major  characteristics  listed  above.  When 
this  has  been  done,  and  done  impartially  and  completely,  we 
will  have  available  for  the  first  time  the  facts  that  are 
necessary  for  the  intelligent  selection  of  an  estimator.  And 
the  selection  should  be  made  in  consideration  of  the  use  that 
will  be  made  of  the  estimate  and  the  characteristics  of  the 
estimator  that  are  important  for  that  application. 

In  many  applications  of  passive  tracking  the  entire  TMA  is 
needed.  For  this  reason  we  believe  that  the  analysis  should  be 
extended  to  TMA  estimators,  as  well  as  range  estimators. 

A  complete  analysis  of  the  type  described  above  seldom  has  been 
attempted,  and  perhaps  never  has  been  done  satisfactorily.  An 
analysis  of  several  range  estimators  (the  Ekelund  and  a  new 
three  bearing-rate  estimator  that  removes  the  bias  in  the  Ekelund 
estimator)  was  made  in  reference  (b) .  Reference  (b)  also  includes 
analyses  of  several  TMA  estimators  (CHURN  and  a  new  one  developed 
by  the  writer),  using  acquisition  probability  as  the  payoff 
function.  While  these  analyses  cover  all  the  major  characteristics 
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of  the  estimators,  they  are  not  extensive  enough  to  be 
considered  complete. 
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